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Abstract 

We discuss a couple of simple quasistatic electromagnetic systems in which the density of elec- 
tromagnetic linear momentum can be easily computed. The examples are also used to illustrate 
how the total electromagnetic linear momentum, which may also be calculated by using the vector 
potential, can be understood as a consequence of the violation of the action-reaction principle, be- 
cause a non-null external force is required to maintain constant the mechanical linear momentum. 
We show how one can avoid the divergence in the interaction linear electromagnetic momentum of 
a system composed by an idealization often used in textbooks (an infinite straight current) and a 
point charge. 

PACS numbers: PACS: 03.50.De 
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I. INTRODUCTION 
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> vS [l|, 3], as well as in a number of articles P, | 
18]. However, the electromagnetic linear momentum 
app earin g in the same kind of systems has comparatively attracted less attention 
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, 0, Q |3| , although it may be even more striking the first time a student is faced with 
it. Furthermore, with the exception of the rather trivial distribution in the example by 
Casserberg j^cj and the mathematically complex calculation in fl7|. in those articles only 
the total electromagnetic linear momentum is computed by using its evolution equation or 
the vector potential. The first goal of this work is to provide some simple examples that 
can be used in the classroom to compute the density of electromagnetic linear momentum, 
as described by the Poynting vector. In our first two examples, the total electromagnetic 
momentum can be computed both by integrating that density and directly by using the 
vector potential. 

In Poynting's theorem one is forced to assume the electromagnetic field carries linear 
momentum in order to get the right balance equation for the total linear momentum. This 
can be illustrated in the particular but simplified context of our examples. In all of them 
the internal electromagnetic forces will not satisfy the action-reaction principle, so that to 
maintain constant the total mechanical momentum one has to apply a non-null external 
force on the system. Then, the only way to keep the balance equation (the total external 
force equals the derivative of the total linear momentum) is to attribute a linear momentum 
to the electromagnetic field. In our examples this balance can be simply checked and one 
may see how the electromagnetic momentum is stored in the field while the system evolves 
starting from a configuration in which some of its elements are separated at infinite distance 
and thus with no electromagnetic momentum. All this may help students to gain confidence 
on their understanding of the meaning of the Poynting vector. 

In Section|n]we review some properties of the kind of quasistatic electromagnetic systems 
we will consider. Then, in Sections 11111 and I1VI we discuss two simple systems in which one 
can easily compute the Poynting vector due to the electric field of a point charge moving 
at constant small velocity and a stationary magnetic field. We explicitly check the balance 
equation of the total linear momentum and also that it may be computed by using either 
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the vector potential or the Poynting vector. We show that the latter method, although more 
difficult in general, has the advantage of showing the spatial distribution of momentum den- 
sity; in fact, it uncovers a curious fact about that distribution in the example of Section HV1 
The total electromagnetic linear momentum of the first model for the system studied in that 
section is infinite. This divergence illustrates the limitations of an idealization often used in 
textbooks (an infinite straight current) and can be avoided in several ways, one of which is 
presented in the same section. In Section IS we analyze another model which can be used 
to illustrate the relationship between the violation of the action-reaction principle and the 
origin of the electromagnetic momentum, as well as the origin of the divergence in the first 
model of Section I1VI We also discuss the energy balance in all the examples. 



II. STATIONARY CURRENT AND POINT CHARGE 

Let us consider an electromagnetic system formed by a static magnetic field B/ created 
by a stationary current / as well as the electric field E g and the magnetic field B g created 
by a point charge q moving with a small constant velocity (which may be null). In the total 
electromagnetic energy there will be a constant magnetic energy associated to the density 
B]/2fi . On the other hand, since the point charge, its fields and the self-contributions to 
the energy density (e £^/2-t-.Bg/2yU ) and to the linear momentum density (e E y x B,j) move 
as a whole, the associated electromagnetic energy and linear momentum are also constant 
(and infinite, so that they must be dealt with by means of some appropriate renormalization 
procedure). 

But from now on we will forget about all these constant self-contributions and consider 
only the interaction energy and linear momentum densities given by the cross terms B^ ■ 
B///i and e E 9 x B/, respectively. 

The total interaction energy is then 

U = qi q -A(r q ), (1) 

where the vector potential in the Coulomb gauge for B; = V x A is evaluated at the charge 
position r q . 

Under the assumption of small constant velocity, the total interaction linear momentum 
can be computed by using the result proved by Calkin [19| (see also 25]) for more general 
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quasistatic systems, which reduces in our case to 



P = gA(r g ). (2) 

Calkin also proves in the more general context that the derivative of the electromagnetic 
linear momentum equals the external force applied on the system to keep constant the 
mechanical linear momentum. One of the goals of this work is to illustrate this kind of 
result ,a S wen as the geeere, _ um baW Q, with som e Partlcularly e asy examples; 
but we will compute not only the total linear momentum but also its density, to gain some 
insight on its spatial distribution. 



III. TOROIDAL SOLENOID AND POINT CHARGE 

Our first example is the system depicted in Figure ^ a point charge is moving with a 
small constant velocity along the symmetry axis of a fixed toroidal solenoid. Since our goal 
is to illustrate in an example a general result, we will make a twofold approximation which 
will allow discussing the physics with a minimal mathematical apparatus. First of all, we 
will assume that the torus is very thin so that the area of the cross section is 5* <C R 2 , 
R being solenoid's main radius. With this assumption the magnetic field created by the 
current I is confined inside the solenoid where its value, in the usual cylindrical coordinates, 
is 

Br = (fi, 3 

2nR K ' 

assuming the solenoid has iV turns. Furthermore, we will assume the point charge located at 
r q = a k moves slowly, so that we will keep in all the calculations only the lowest order term 
in an expansion in powers of (3 = a/c. In particular, we will check below that this allows 
ignoring the change of the current intensity on the solenoid due to the changing magnetic 
flux created by the charge, because the induced electromotive force is proportional to f3 2 . 
Moreover, this approximation reduces the exact expressions for its electric and magnetic 
fields at a generic point of coordinates (p, ip,z), 

E = (1 - /3 2 ) g pp+(z-a)k 

11 ' 47Te [(l_/32)p2 + ( z _ a )2]3/2' 

B q = ^Tg X Eg, (5) 
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to the easier 

E, = -^- P * +( '~ a) 3 k /a , (6) 

47Te [p2 + ( z _ a ) 2] 3/2 
9 4VT [p2 +( ^_ a) 2]3/2^ 

A. Energy balance 

The total flux of the changing magnetic field B g across the solenoid located at 
(p = R, z = 0) is 



poqNRSd 
4tt (i? 2 + a 

so that the induced electromotive force is 



/£= (10) 



rf$ _ . d$ 3p qNRSaa 2 _ 3qNRSa a 2 
' ~ltt ~ ~ a l^ ~ 4vr (R 2 + a 2 f 2 ~ 4vre (R 2 + a 2 f 2 c 2 ' 

The dependence on /3 2 of the last expression explicitly shows that, in our approximation, 
we could ignore induction when computing the lowest order of the magnetic field created by 
the solenoid. The power given to the solenoid is, then, 

3p qNIRSaa 2 
4tt (R 2 + a 2 

Since (in our low-velocity approximation) the solenoid creates no magnetic field outside 
the torus, the point charge feels no magnetic force: no external force is necessary to keep 
its constant velocity. But, due to the magnetic field created by the charge, the solenoid will 
feel a magnetic force, so that an external force is necessary to keep it at rest. Since it is 
applied on a fixed solenoid, the external force performs no work. Furthermore, the kinetic 
energy is constant. So, where does power (fTU|) come from? 

The answer is, of course, that the electromagnetic energy is not constant: it has a time- 
dependent component due to the mixed electromagnetic field, which can be computed by 
using (^) or directly from the density of electromagnetic energy: 

1 f TD ~D J~\ r RS I TD TD J.. P qNIRSd 



U = qr- A(r,) = — / B q -B l( iV= I Bj-B ff cfy> 

Un JT Un JT 
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PqJt - p Jt ' 4n{R 2 + a 2 f 2 ' 

We can now see that the derivative U is just the opposite of the power given to the solenoid 
by the changing magnetic field B g , so that the conservation of the total energy holds: 

Ie + U = 0. (12) 
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B. Linear momentum balance 



U=j>dm-B q = ^ (15) 



As mentioned above, the point charge feels no magnetic force and there is no external 
force applied on it. However, we can easily check that the action- reaction principle is not 
satisfied, for there is a magnetic force on the solenoid due to B g . 

Let us consider the current turns located between the angles p and p> + dip. We may 
think about them as an equivalent elementary turn of current 

dI = NI% (13) 

so that its magnetic dipole moment is 

dm — dl S (p — dp (p. (14) 

Incidentally, we can now check that the electromagnetic energy (fTTj) can also be computed 
by summing up the energies of these elementary dipoles: 

fj, qNIRSd 
lc H 4vr (R 2 + a 2 

On the other hand, as a consequence of the magnetic field created by the moving charge, 
there is a force acting on the turn, which can be easily calculated: 

fipqNISd p 

dmB '- ar> v + iz-ayr*- (16) 

V qJ\ p=R)Z=0 (i?2 +Q 2)5/2 V V ) 

The total magnetic force on the solenoid is now readily computed: 

3fi qNISRad 
lc 4vr (R 2 + a 2 

Let us now assume that the particle came from infinity with a small constant velocity 
a < 0. The particle feels no force, but to keep at rest the solenoid one has to apply on it an 
external force F ext = — F. The total impulse of this external force is readily computed: 

da _ f°° 1 ^ ^ _ fi qNIRS 
4tt (R 2 + a 2 

At first sight one could think it is contradictory to have a net external impulse while the 
mechanical linear momentum remains constant along the whole process. To solve this puzzle 



F=f_dF= 5/2 k. (18) 



/ P_ dt = f P_ $ = f i F da = 7f 7 ^ k. (19) 
7 Joo a J a a 4tt R 2 + a 2 ) 3 ' 2 



and keep the general result that states the total impulse applied on a system equals the change 
of its total linear momentum, one has to take into account the linear momentum stored in 
the electromagnetic field. 

The interaction electromagnetic linear momentum density in this setup is confined inside 
the solenoid, where its value is 

r - e F v R - ^° qNI a P+ Rk / 9n \ 
G-6 ^xBj- Q ^ 2p ^ (20) 



8n 2 R (R2 



a 



We see that the linear momentum density has two components which vanish at infinity. 
The one along the symmetry axis OZ is always positive and decreases monotonically from 
a = to a — > ±00. The radial component has the sign of a and reaches its maximum 
absolute value for a = ±R/y/2. By symmetry this radial component will give a total null 
contribution to the linear momentum. In fact, the latter is easily computed, because in our 
approximation the integral extended to the inside of the torus T reduces to an integral along 
its central line C: 

P= / GdV = SR<( G#= mNIRS ^ 2 k. (21) 
Jt Jc r 4vr {R 2 + a 2 f 2 

The physical picture is now clear: the electromagnetic linear momentum ({21} stored 

inside the solenoid is a direct consequence of the impulse of the external force F ext = — F 

applied on the solenoid to keep it at rest. If the particle came from infinity (where the 

electromagnetic linear momentum was null) with a small constant velocity a < 0, the total 

impulse of the external force equals the total electromagnetic momentum, because from (JT5J) 

and (J21|) we have 

jF ext dt = P (22) 

or, equivalently, F ext dt = P. 

On the other hand, as a consequence of (|14|). the vector potential created at charge's 

location by the elementary turn (p, tp + dip) is 

^0 dm x (r q - r) fx NIS a p + Rk 

dA (r„) = ■ — = — — ^ dip. (23) 

4vr | r , - if 8tt 2 (#2 + a 2)3/2 

It is remarkable that, from the point of view of result ©, the linear momentum inside 
the elementary turn may be understood as created only by that turn, independently of the 
remaining turns in the solenoid: 

dP = GdV = GSRdp = qdA(r q ). (24) 



Of course the same is not true outside the torus, where the linear momentum density is zero 
due to the sum of contributions by all turns. 

We can now check directly that (J2J) gives ([21)1. because from the we find the total 
vector potential: 

This direct method to compute the linear momentum is not easier than the calculation 
involving the Poynting vector — unless (|25|) is computed (by checking it gives the corretct 
B/, for instance) without calculating (J23)) — , which would give us no hint on the presence 
of a radial density of linear momentum. 



C. Angular momentum balance 

Since the dipole moment ()14|) is parallel to the magnetic field created by the charge, there 
is no torque on the turn, dm x B, = 0, which is consistent with the constant (null) value 
of the angular momentum. In fact, the mechanical angular momentum — with respect to 
the center of coordinates — is null and the density of electromagnetic angular momentum is 
parallel to the solenoid, 

HoqNIR A 

rxG = Tj^V, (26) 

8tt 2 {B 2 + a 2 ) 3/2 

so that there is no electromagnetic angular momentum: 

L = jf r x GdV = Sr£ t x G d V = 0. (27) 

IV. STRAIGHT CURRENT AND POINT CHARGE 

Let us now consider a fixed infinite straight wire conducting a constant current / and 
a point charge directly moving away (or towards) the wire along a perpendicular direction 
with constant small velocity. If we choose cylindrical coordinates around the OZ axis of 
Figure El the charge position is given by 



r 



a i = a (cos cp p — sin tp (p) . (28) 



In the following we will neglect the wire radius and, as in the previous section, we will only 
keep at each step the lowest order in (3 — d/c, so that, at a generic point r = p p + zk, the 
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electric field of the charge and the magnetic field created by the current are, respectively, 



q (p — a cos (p)p + asmtpfi + zk 
47re (p2 _|_ a 2 _|_ z 2 _ 2 pa cos 



A. External force 

By symmetry the magnetic field B g created by the charge exerts no net force on the wire, 
but, again, the action- reaction principle is violated because there is a magnetic force exerted 
on the point charge by the magnetic field B/, so that to keep constant the mechanical linear 
momentum one has to apply on the point charge an external force opposite to the magnetic 
one: 

F ext = -qr q x B 7 . (31) 

Thus, we need a changing electromagnetic linear momentum whose time derivative equals 
the total external force (|31|) applied on the system. 

B. Linear electromagnetic momentum 

The cross term in the linear momentum density is easily computed from (|29~j) - (|3()|) : 



G = e E q x Bj = M£ -^+(p-acos^)k 

8tt 2 p ( p 2 + a 2 + z 2_ 2pa cos <^) 3/2 



Since dr oc G along the current lines of the linear momentum density, their differential 
equations read as follows: 

^= dZ , &p=0. (33) 
— z p — acosip 

The first equation is readily integrated to give (p — acos<^) 2 + z 2 = C 2 , so that the current 
lines of this density field are arcs of circumference with center at a cos ip p lying in planes 
tp = const. Some of them are depicted in Figure H3 By using a well known identity, we get 

V • G = e B 7 • (V x Eg) — e E 9 ■ ( V x Bj) = -^E q ■ j, (34) 

where j is the current density and we have used V x E, = and VxB; = p j- So, as can 
also be proved directly from (j32jl . the divergence of G vanishes except inside the conducting 
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wire, where it is positive for z < and negative for z > 0. This in turn explains why its 
current lines are closed or go from a wire point with z < to another wire point with z > 0. 

Due to symmetry, the radial component of density (|32j) will vanish when integrated with 
respect to z, so that the total electromagnetic linear momentum due to the interaction of 
the electric field of the point charge and the magnetic field created by the current is 

td In a a a 1 f°° a f 277 a p a, cos (p p ql r°°dp 

P = / Gpdpdpdz = — — k / dp dtp — = — — k / — . 35 

J Air z Jo Jo p + a — 2pa cos (p 2n J a p 

We have used here the result (jAlj) . which shows a curious fact: there is no net contribution 
to the electromagnetic linear momentum for p < a, i.e., inside the cylinder centered at the 
current which goes through the charge. Although it is clear from (J32|) — and from Figure 0J— 
that this is only possible for p < a, because only there does the sign of G z change when p 
varies, one might wonder whether this exact cancellation happens only in the low velocity 
approximation we are using. The answer is that the same cancellation occurs when the 
exact electric field of the point charge is used. This can be easily checked, because the linear 
momentum density in the relativistic case is 

/i g/(l-/? 2 ) -z p+(p- a cos p)k 

— ^r n — — — : — - — z~ — : z : . ~„ ^ w 



in*p 



p 2 + a 2 + (1 — (3 2 ) z 2 — 2pa cos <p — (3 2 p 2 sin <p 



so that using the integral (jA2|) we get 



P = / G p dpdpdz = — — k f dp [ dp 1 r 



An 2 Jo Jo p 2 + a 2 — 2pa cos ip — (3 2 p 2 sin <p 

M/ k f^. (37) 



2lT Ja p 

On the other hand, the contribution to the momentum (J35|) for p > a has a logarithmic 
divergence, as one could expect from the slow decrease of the magnetic field ()3())1 at infinity. 
This is just a consequence of the well known limitation of the model we are using: an infinite 
straight current is a good approximation to compute fields near the wire, but it may lead 
to meaningless results for far fields, which enter along with near fields in the total value of 
the electromagnetic momentum. This limitation can be easily avoided if we consider that 
the return current goes along a coaxial conducting cylinder of radius p = R > a, so that 
it encloses both the wire and the particle. Then B/ and the Poynting vector will vanish 
outside the cylinder and instead of fl35(l we have the finite 

« m/ k f^ = ^ln^k. (38) 
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The derivative of P is well defined in both models (and directly related to the change 
due to the motion of the surface p = a with separates the regions with null and non-null 
contribution) : 

P = a— = -^1— k = -qv q x Bj = F cxt . (39) 

This is the right evolution equation for the linear momentum: the net external force applied 
on the system equals the time derivative of the total linear momentum. 

The vector potential of the magnetic field created by the currents through the wire and 
the cylinder is in our approximation 

A=^9(R-p)\n*k, (40) 

Z7T p 

where the Heaviside function is defined as usual: 

1, if x > 0: 



6{x) = { 



(41) 
0, if x < 0. 



So, if we were interested only in the total electromagnetic linear momentum we could have 
used (I4U|) for p = a in (J2J) to compute (|38|) more easily, but then we would have known 
nothing about the curious distribution of its density. 

In this example the magnetic force applied on the charge performs no work, so both the 
mechanical and the electromagnetic energies are constant. In fact, it is easy to see that the 
density of electromagnetic energy l///oB 9 -Bj is antisymmetric with respect to the reflection 
z <-> —z and, in consequence, the total electromagnetic energy due to the interaction of B g 
and B/ is zero, which is also a direct consequence of ([TJ because the charge velocity and the 
vector potential are perpendicular. 

By taking the limit R — > oo in (|3*Hj) we recover the divergence in ()35|) . but the latter can 
also be understood by considering another family of models that recovers in the appropriate 
limit the one we have used at the beginning of this section. 



V. CIRCUIT AND POINT CHARGE 

In this section we will consider the circuit of Figure 0] and the point charge moving with 
constant velocity along a symmetry axis. If we take the double limit L, d — >• oo we recover 
the system of Section HVl By symmetry, the flux of the magnetic field created by the charge 
across the circuit is null, so that there is no induction. 
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In this case also the action-reaction principle is violated, for the external forces that have 
to be applied on the charge (to keep constant its velocity) and on the circuit (to keep it at 
rest) are not opposite to each other. By symmetry, the magnetic force exerted on circuit 
sides 1 and 3 vanish, while those applied on sides 2 and 4 are equal, so that the total magnetic 
force on the circuit is 

^oqld I a + d a 



2tt 



(42) 



L^J(a + d) 2 + L 2 Wa 2 + L 2 / 
as one can easily check. It is not difficult to compute the magnetic field created by each 
circuit side on the point charge. The force exerted by the magnetic fields of arms 2 and 4 
is just — Fj, but there are also the forces due to the magnetic fields created by 1 and 4: the 
total magnetic force on the charge happens to be 

fioqld ( L L a a + d 



2tt y a Va 2 + L 2 ( a + d)^J( a + d) 2 + L 2 Wa 2 + L 2 L^/(a + d) 2 + L 2 

'(43) 

In consequence, to keep at rest the circuit and constant the charge velocity, external forces 
must be applied on both elements, so that although the mechanical linear momentum is 
constant, the net external force on the system does not vanish: 

F ext = - (F, + F a ) = ( - L , L ) k. (44) 

K 1 q> 2vr [aVa^+I? (a + d)^(a + d) 2 + L 2 ) 1 1 

Again we need a changing electromagnetic linear momentum to have the right evolution for 
the total linear momentum. 

The electromagnetic linear momentum due to the interaction of fields created by circuit 
and charge can be easily computed by using (J2J) and 



Hoi / dr 



A(r ' )= ^i^ (45) 

where C is the circuit. By symmetry, the contribution from sides 2 and 4 to (J45)) are opposite: 
A-2 (r g ) + A 4 (r g ) = 0. The contribution from side 3 is 

A 3 (r 9 ) = -^L k [ L , dz = arcsinh — L — k. (46) 

47T J-L /( a + d) 2 + Z 2 27T \a + d\ 1 > 



Since the contribution from the remaining side is Ai = —A3 , the vector potential (|43|) 

d — s-O 

at the location of the point charge is 

A (r„) = (arcsinh — — arcsinh -) k. (47) 

v ql 2tt \ \a\ \a + d\J v ' 
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Now we can check directly the evolution equation for the total linear momentum: 

P = Jt [qA (Tq)] = q ^ (r?) = Fext - (48) 

On the other hand, from result ()47|) it is easy to understand the origin of the logarithmic 
divergence of the vector potential and the linear momentum qA (r q ) in the limit case of the 
beginning of Section HV| because if we compute in any order the double limit we get: 

a + d 



lim A (r g ) = k lim In ^— — — — = k lim In 

d^oo, L~ *oo 2,71 L^oo fl 2ll d-^oo 



(49) 



As in Section II VI the external forces perform no work and the charge velocity and the 
vector potential are perpendicular, so that the total electromagnetic energy due to the 
interaction of B g and B/ is zero. 

Acknowledgments 

This work was supported by The University of the Basque Country (Research 
Grant 9/UPV00172.310-14456/2002). 

APPENDIX A: A COUPLE OF USEFUL INTEGRALS 
To compute (j3*5Jl . we have used 



/ 27r p- acosip , 27T ^ 

^= 2T 2 9 d<p=— 9(p-a) (Al) 

o p + a — 2pa cos (p p 



and in ()37|) one needs the more general 

f 2n p-acos(p , 271-^. 

^= 2^2 O ^22-2 d<P= n m d (p- a )- A2 

Jo p z + a z — lap cos (p — p z p z sin (p pv 1 — p z 
Integral (jAljl — as well as others appearing in the previous pages — can be computed by 



using tables [28j or a computer algebra system, but both methods are far less helpful in 
order to calculate (|A2|) . In fact, the easiest way to prove results (|A1J) and (|A2J) is to consider 
in the complex plane the ellipse C described by 



z = (p cos ip — a) + ipJ 1 — (3 2 sin <p, (0 < (f < 2n) (A3) 
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for constants p > 0, a > and < (3 < 1. Then, since the pole z = lies inside C for p > a, 
we have 

h = \mj c ^= p^\-^h = 2n9(p - a), (A4) 

which is equivalent to ()A2j) and reduces to (jAljl for (3 = 0. 

For students who do not know about complex integration, one might use an equivalent 
but more pictorial method. By eliminating the parameter < <p < 2tt, it is easy to check 
that the equations 



x = pcos(p — a, y = p\J 1 — (3 2 simp (A5) 

describe an ellipse with center at (—a, 0), horizontal major semi-axis of length p and eccen- 
tricity P, as depicted in Figure El If we consider the polar coordinates (r, 6) of an ellipse 
point P, the integral (jA4|) is just the full change of the polar angle 6 in a complete counter- 
clockwise turn of P around the ellipse: 



A# = f d6 = f darctan— = f (iarctan = pJl — (3 2 I 2 . (A6) 

Jc Jc x Jc pcosLp — a v 

It is not difficult to convince oneself that A8 equals 2ir if the origin is inside the ellipse (i.e., 
if p > a) and otherwise. 
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FIGURE CAPTIONS 



Fig 1 Toroidal solenoid and point charge 

Fig 2 Straight conducting wire and point charge 

Fig 3 Current lines of the linear momentum density (J32|) . In the upper left corner the lines 

3a 71 71 37T 

passing through (p, <p, z) = (— , <p, 0) for ip = 0, — , — , — . In the other three corners 

TT 7T 

some lines in the planes p = 0, — , — are drawn, including — in bolder line — the one 
also displayed in the upper left corner. The intersection between the plane and the 
cylinder p = a appears as a couple of broken vertical lines. 

Fig 4 Circuit and point charge 

Fig 5 Ellipse with center at (—a, 0) and semi-axes equal to p and py/1 — (3 2 
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FIG. 1: Toroidal solenoid and point charge. 
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AguirregabiriaFigl 
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. 2: Straight conducting wire and point charge. 
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AguirregabiriaFig2 




FIG. 3: Current lines of the linear momentum density (|H2j) . In the upper left corner the lines 
passing through (p,(p,z) = (— ,tp,0) for cp = 0, — , — , — . In the other three corners some lines 

7T 7T 

in the planes = 0, — , — are drawn, including — in bolder line — the one also displayed in the 

upper left corner. The intersection between the plane and the cylinder p = a appears as a couple 

on 

of broken vertical lines. 
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FIG. 4: Circuit and point charge. 
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AguirregabiriaFig4 
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FIG. 5: Ellipse with center at (—a, 0) and semi-axes equal to p and p^J\ — (5 2 . 
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